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We consider open many-body systems governed by a time-dependent quantum master equation 
with short-range interactions. With a generalized Lieb-Robinson bound, we show that the evolution 
in this very generic framework is quasi-local, i.e., the evolution of observables can be approximated by 
implementing the dynamics only in a vicinity of the observables' support. The precision increases 
exponentially with the diameter of the considered subsystem. Hence, the time-evolution can be 
simulated on classical computers with a cost that is independent of the system size. Providing error 
bounds for Trotter decompositions, we conclude that the simulation on a quantum computer is 
additionally efficient in time. For experiments and simulations, our result can be used to rigorously 
bound finite-size effects. 
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I. INTRODUCTION 

In Lorentz-invariant theories, a maximum speed for the 
propagation of information is, by construction, the speed 
of light. In nonrelativistic quantum theory, the existence 
of a maximum propagation speed results more indirectly 
and for different reasons. For nonpathological models, 
this maximum speed is much smaller than the speed of 
light. The seminal paper by Lieb and Robinson [1] and 
further contributions like [2-13] cover isolated systems. 

Here, we consider the evolution of a more general 
and, experimentally, extremely relevant class of systems - 
open quantum many-body systems governed by a quan- 
tum master equation [14, 15] with short-range Liouvil- 
lians that are allowed to be time-dependent. Promi- 
nent experimental examples are presented in Rcfs. [16- 
20], and recent theoretical advances on quantum com- 
putation, noncquilibrium steady states, and phase tran- 
sitions in open systems can, for example, be found in 
Refs. [21-24]. Going beyond the existence of a finite 
maximum propagation speed and the existence of a well- 
defined thermodynamic limit [1, 25], we show that the 
time-evolution of such systems is quasi-local. This means 
that, up to an exponentially small error, the diameter 
of the support of any evolved local observable grows at 
most linearly in time, or, put differently that the evo- 
lution of the local observable can be approximated to 
arbitrary precision by applying the propagator of a spa- 
tially truncated version of the Liouvillian; Fig. lb. For 
the special case of isolated systems, where the evolution 
is given by a unitary transformation, the corresponding 
question has been addressed in Ref. [9]. As a tool for 
the proof of quasi-locality, we derive and employ a Lieb- 
Robinson-type bound very similar to the recent results 
of Poulin [26] and Nachtergaele et al. [25] . All constants 
in the bounds are given explicitly in terms of the system 
parameters. 

The quasi-locality of Markovian quantum dynamics 
has several crucial consequences. It implies that the evo- 
lution of observables with a finite spatial support can be 



simulated efficiently on classical computers, in the sense 
that the computation cost is independent of the system 
size, irrespective of the desired accuracy. This can for ex- 
ample be exploited in an exact diagonalization approach 
for a sufficiently large vicinity of the support of the con- 
sidered observable; Fig. lb. For more sophisticated sim- 
ulation techniques, we provide, in extension of Ref. [27], 
error bounds for Trotter decompositions [28] of the sub- 
system propagator into a circuit of local channels; see 
Fig. lc. The Trotter error is polynomial in the time, at 
most linear in the size of the time step, and can hence be 
made arbitrarily small. Importantly, the subsystem Trot- 
ter decompositions allow for the efficient simulation of the 
time-evolution on a quantum computer as envisaged by 
Feynman. For any required accuracy, the simulation can 
be implemented with a cost that is independent of the 
system size and polynomial in the time. 

Experimental and numerical physicists who study 
nonequilibrium systems can use our result on quasi- 
locality to rigorously bound finite-size effects. This is for 
example relevant for experiments with ultracold atoms 
in optical lattices [29] and numerical investigations em- 
ploying time-dependent density-matrix renormalization 
group methods [30-33]. 

II. SETTING 
A. Lattice and equations of motion 

Let us consider lattice systems, where each site z e A 
is associated with a local Hilbert space rl z . Subsystem 
Hilbert spaces are denoted by 

H v ■= (g) H z Vy cA and H ■= H\. 

z£V 

Let p{t) denote the system state at time t. Markovian 
dynamics of an open quantum system, i.e., the evolu- 
tion under a linear differential equation that generates a 
completely positive and trace-preserving map for p, can 



2 



always be written in the form of a Lindblad equation [34- 
36] 

d t p = -i[H, p]+J2 UupLl - \{LlL vP + P LlL v )\ , 

where the arbitrary Lindblad operators L v and the Her- 
mitian Hamiltonian H may depend on time. This equa- 
tion captures, for example in the framework of the Born- 
Markov approximation, the evolution of a system that 
interacts with an environment [14, 15] and isolated sys- 
tems as a special case. Let us switch from the Schrodingcr 
picture, where expectation values are evaluated according 
to (0} s ^t = Tr[p(t)0] with p(s) — p, to the Heisenberg 
picture, where (0) s _> t = Tr[pO(s)] with 0(t) = O. The 
corresponding timc-dcpcndcncc of an observable O(s) £ 
B(H) is then given by the quantum master equation 

d s O(s) = -£(s)0(s), 

where £(t) e B(B(H)) is a super-operator, the so-called 
Liouvillian, with the Lindblad representation 

CO = i[H,0] + UlOK - \{L\L v O + OLtL v )\ . 

The set of Liouvillians with spatial support V C A will 
be denoted by L y c B(B{U V ))- 

B. Short-range Liouvillian 

In order to be able to use Lieb-Robinson bound tech- 
niques, we need to restrict ourselves to Liouvillians with 
norm-bounded short-range interaction terms. Let us 
hence assume that £ is a sum of local Liouville terms 
lz with norm bound maximum range a, and a max- 
imum number Z of nearest neighbors [47]. Specifically, 



£(t) 


= M*)> £ z (t) e L z , 


(1) 




ZcA 




\t\ 


■= sup \\l z (t)\\, 


(2) 




t,ZcA 




a 


:= sup diam(Z), 


(3) 




Z:£ z ^0 




Z 


:= max \{Z' C A \l z > i- 0, Z' n Z ^ 0}|, 


(4) 



where diam(Z) := maxx^gz d(x, y) is the diameter of 
Z and d is a metric on the lattice A. In Eq. (2), we 
have used the super-operator norm defined by ||T|| := 
sup 06B ( W ) \\TO\\ 1 1|0||. In the Heisenberg picture, this is 
the physically relevant norm as induced by the operator 
norm ||0||; see Appx. A. For notational convenience, 
we define for every subsystem V C A the corresponding 
extension V, volume Vol(V), and truncated Liouvillian 
Cy, 

V:=\J z .. iz ^ Z, (5) 

Vol(V):=\{ZcV\£ z ^0}\, (6) 
£v(t):=E Zc vWt)- (7) 



C. Propagators 

Propagators ry(s,t) are super-operators that map ob- 
servables to time-evolved observables. They are defined 
as the unique solutions of 

d s T V (s,t) = -£ v (s)T V (s,t), T V (t,t) = id V s < t . (8) 

With t(s, t) := ta(s, t) one has indeed O(s) — t(s, t)0{t). 
Propagators obey the composition rule r(r, s)r(s,t) — 
r(r,t) y r <s<t- As discussed in Appx. B, the derivative 
with respect to the second time argument is given by 

d t T V (s,t)=T V (s,t)£ v (t), (9) 

and the propagators are norm-decreasing, 

\\r(s,t)0\\ <\\0\\ \f£eh A ,s<t,OeB(n). (10) 

III. QUASI-LOCALITY OF THE EVOLUTION 

Given an operator Oy G Billy) with support Y C A, 
we would like to show that the exactly time-evolved op- 
erator r(r,t)Oy with r <t can be approximated by the 
evolution with respect to a spatially truncated Liouvil- 
lian, i.e., by Ty(r,t)Oy with Y C V C A. Indeed, our 
main result, Theorem 2, states that the approximation 
error is exponentially small, in the distance of A \ V to 
the time-r slice of a space-time cone originating from the 
operator's support Y at time t; see Fig. lb. More pre- 
cisely, the error decays exponentially in d(Y, A \ V) / 'a — 
v ■ (t — r), where d(X, Y) := mf xe x,yeY d(x, y) is the dis- 
tance of two subsystems X, Y C A, and v = eZ\£\ is the 
so-called Lieb-Robinson velocity. 

To prove this, we can write the difference of the evolved 
operators in the form 

r(r, t)Oy - Ty{r,t)Oy 

= - J &sd s [T V (r,s)T{s,t)]Oy 

= [ AsT v {r,s)[£{s) - £ v {s)]T{s,t)Oy 

Jr s v ' 

=£a\vO) 

due to the fundamental theorem of calculus and Eqs. (8) 
and (9). Using the triangle inequality and the fact that 
the propagators are norm-decreasing, it follows that 

\\r(r,t)0 Y - T v (r,t)Oy\\ 

< f ds ¥x{ S )r{s,t)Oy\\. (11) 

In the case of unitary dynamics (ix(s)O = i[h x ,0]), the 
integrand would be of the form || [hx , t(s, t)Oy] || , and the 
standard Lieb-Robinson bound [1-5] would be applicable. 
To proceed in our more general case, however, we use a 
Lieb-Robinson bound for Markovian quantum dynamics, 
similar to recent results in Refs. [25, 26]. 
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Figure 1: (a) An evolved local operator r(s,i)Oy behaves almost like the identity outside its associated space-time cone, (b) 
Approximating r (s, t)Oy by application of subsystem propagators to Oy ■ The errors decrease exponentially with the subsystem 
sizes, (c) For one-dimensional systems, approximating r(s, t)Oy by a Trotter decomposition yields an error scaling as (t — s) 2 At. 
Note that the Trotter circuit can be trimmed off at the boundary of the Lieb-Robinson space-time cone. 



Theorem 1 (Lieb-Robinson bound for Markovian 
quantum dynamics) 

Let the Liouvillian C(t) — X^zcA^-^(^) f or the lattice A 
be of finite range a, with a finite maximum number Z of 
nearest neighbors, and \t\ as defined in Sect. II. Also, let 
K x £ h x , Y £ B(U y ), and r <(eR. Then 

\\JC x r{r,t)0 Y \\ < Vx.y \\Oy\\ e^-^-'W", 

(12) 

where v '■= eZ\£\ and Vx,y '■— min{ Voi ^ x ^ ; v °K y ) } . 

The proof is given in Sect. V. The theorem tells us 
that an evolved observable r{r,t)Oy remains basically 
unchanged when we evolve it with respect to a Liouvil- 
lian that is supported at a distance R^> v ■ (t — r) away 
from Y, i.e., that T{r,t)Oy behaves like the identity out- 
side the corresponding space-time cone. In the special 
case JCxO = i[Ox,0], Eq. (12) yields a Lieb-Robinson 
bound for \\[O x ,T(r,t)0 Y ]\\ as in Ref. [26]. 

This theorem can now be employed to proceed from 
Eq. (11) in our proof of quasi-locality. Let us restrict 
ourselves to the typical case of Liouvillians C(t) for which 
the number of terms £x(t) with distance d(y,X)/a £ 
[n, n+ 1) from any site y £ A is bounded by a power law, 

\Rn, y \ < Mn K V ye A : „ eN+ , (13) 

R n>v :={XCA|I^0, e [n,n+l)}, 

for some constants M, k > 0. Now, choose a point yo £ Y 
that is closest to A \ V, i.e., d(y , A \ V) = d(Y, A \ V). 
With D := \d(Y, A\V)/a'}, we can exploit that the sup- 
port of every term in Ca\v 1S element of exactly one of 
the sets R n ,y with n > D, to obtain 

\\T(r,t)0 Y -T v (r,t)0 Y \\ 

^ E E / ds\\lx(s)r( S ,t)0 Y \\ 

oo 

< Mn K \£\ \\0 Y \\ / Ase"- {t - r) - n 

n—D 

p v-(t—r) 00 

< M\£\ \\O y \\ V n K e~ n . 

v ^ — ' 

n=D 



In the second step, Theorem 1 and Vxy < Vol(X)/Z < 
1 have been used. With the bound Y^=d nKe ~ n < 
2eD K e~ D Vd>2 K +i from Appx. D, we arrive at the cen- 
tral result of this work: 

Theorem 2 ( Quasi-locality of Markovian quan- 
tum dynamics) 

Let the Liouvillian C(t) = X)zcA^ z W f or ^he lattice A 
be of finite range a, with a finite maximum number Z 
of nearest neighbors, and \£\ as defined in Sect. II. Fur- 
ther, let constraint Eq. (13) be fulfilled for some constants 
M,K > 0. Also, let Y C V C A, Y £ B{U Y ), and 
r <t£R. Then one has with D := \d(Y, A \ V)/d\ 

\\r(r,t)0 Y - Ty{r,t)0 Y \\ 

< 2M || 0y || DKeV . {t - r )-D Vd>2k+1 , (14) 

where v is the Lieb-Robinson speed from Eq. (12). 

The full dynamics can be approximated with exponential 
accuracy by subsystem dynamics. In a sense, the con- 
straint Eq. (13) requires the lattice to have a finite spa- 
tial dimension. A 2?-dimensional hypercubic lattice with 
finite-range interactions fulfills Eq. (13) with n = T> — 1. 
An interesting observation is that short-range models on 
a Bethe lattice [37] have a finite Lieb-Robinson speed ac- 
cording to Theorem 1 but do not fulfill Eq. (13) and are 
thus not covered by Theorem 2. For such systems, it is, 
hence, conceivable that a quench of the Liouvillian start- 
ing at time t = with a distance of at least aD from 
some point y causes a perceptible effect at y for a time 
t* < D/v. 

IV. TROTTER DECOMPOSITION OF THE 
EVOLUTION 

The quasi-locality of the dynamics, Theorem 2, im- 
plies that the evolution of obscrvables with a finite spatial 
support can be simulated efficiently on classical comput- 
ers, in the sense that the computation cost is indepen- 
dent of the system size, irrespective of the desired ac- 
curacy. However, exploiting this in an exact diagonaliza- 
tion approach that stores the approximated time-evolved 
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observable T v (r,t)Oy in a full basis of % v exactly, re- 
quires resources that are exponential in the size \V\ of 
the considered subsystem. There are more sophisticated 
numerical techniques, e.g., one can use matrix-product 
operators [38-40] for the representation of (an approxi- 
mation to) ry(r,t)Oy or sampling algorithms. In such 
schemes, it is typically not possible to address the differ- 
ential equation for Ty(r,t)Oy directly, but one can use 
Trotter decompositions [28] instead, where propagators 
ry (r, t) are decomposed into a circuit of local (diamcter- 
a) channels. 

Using the quasi-locality, Theorem 2, and techniques as 
in Ref. [27], we can derive a Trotter decomposition with 
an error that is polynomial in time, at most linear in the 
time step, and, in extension of Rcf. [27], system-size inde- 
pendent. Furthermore, implementing such a Trotter cir- 
cuit on a quantum computer [27] yields a simulation that, 
additionally to being independent of the system size, is 
efficient in time. In this case, the physically relevant 
norm for super-operators T is the subsystem-seminorm 



||T|| y := sup \\TO Y \\/\\Oy 
OyeefHy) 



(15) 



Theorem 3 (Efficient Trotter decomposition of 
time-evolved observables) 

With the preconditions of Theorem 2, a sequence of 
times to < t\ < • • • < t x and a sequence of sub- 
systems Y c Vi C Vi C • • • C Vjv C A such that 
D n := \d(Y, A\ V n )/a] > 2k + 1 V„ ; the Trotter decom- 
position 



N 



f := I! j I T z(tn-l,t n ) 

"=1 ZC?„:I 2 ^ 



(16) 



into propagators tz for local Liouville terms iz approxi- 
mates the full system propagator r(to, ijv) U P t° an error 



N 

\r(t ,t N ) -f\\ Y < J2( 2 -¥D^-^- D -+s n ), 

n=l 

\2cnr^fTr M D\ 2 Jt n -t n -!)\e\ 



s n := (t n - t„-i) 2 ^ol(V n )|*| V*""*"- 1 '"'! (17) 

with the Lieb-Robinson speed v from Eq. (12). 

In the Trotter decomposition f , we used the convention 
n^Li = T\T 2 . . . Tjv, and the ordering of the channels 
tz in the second product of Eq. (16) can be chosen arbi- 
trarily. As in Ref. [27], one can use averaged Liouvillians, 
i.e., Tz(r,t) i y e S r Asl z(s)^ without changing the scal- 
ing of the error bound. Choosing a constant time step, 
t n = nAt, and subsystems V n such that D n = D +vnAt, 
for sufficiently large D , the bound (17) is dominated by 
the Trotter errors s n . The subsystems can be chosen such 
that diamV^ < diam(y)+a_D„; sec Fig. lc. For this case, 
the total error is in O (Ai(diam(Y)/<z + D + vt) K+2 ). 
Higher-order Trotter-Suzuki decompositions [41] can be 
used to further improve the scaling in At. 



To prove Theorem 3, one can first apply Theorem 2, 
the inequality ||TiT 2 - TiT 2 || < ||Ti|| ||T 2 - f 2 || + ||Ti - 
Ti||||!Z2||, and Eq. (10) iteratively N times, to obtain 



N 



|r(t 0) t N ) - r v \\ Y < D«e-^-*°)- D " (18) 



n=l 



with t v := Yln=i T v n (tn-i,t n )- For every time-step 
propagator Ty (t n -i,t n ), we can then employ a Trotter 
decomposition similar to Ref. [27], yielding 

\\ry(r,t)- J] r z (r,t)||y 

zcv,e z ^o 

< {t-r) 2 ZYol(V)\£\ 2 e^ rW . (19) 

See Appx. E for details. Combining Eqs. (18) and (19) 
with the triangle inequality proves Theorem 3. 



V. PROOF OF THEOREM 1 

With an argument similar to those in Refs. [1-5, 25, 
26] , we want to bound the norm of the operator 



G(r) := JC X T{r,t)0 Y 



(20) 



under the preconditions of Theorem 1. G is the solution 
to the final value problem G(t) = JCxOy, 

d r G(r) = -K x C{r)T{r,t)0 Y 

= -£ A \x(r)G(r) - /C x C x (r)r(r,t)Oy, 

due to Eq. (8), jC = £ x + £ A \x[Eq. (5)], and IC X £ A \ X = 
£\\ X JC X for all Liouvillians JC X G h x . As can be 
checked by differentiation, a corresponding integral equa- 
tion for G(r) is 

G(r) - T AXX (r,t)G(t) 

+ / dsT A \ x (r,s)/C x £ x (s)T(s,t)Oy. 

J r 

Using the triangle inequality, the norm-submultiplica- 
tivity, and the fact that the propagators are norm- 
decreasing, this yields the bound 

l|G(r)||<||G(t)|| + ||^|| [* ds \\£ x (s)T(s,t)0 Y \\ 



< \\G(t)\\ + \\IC X \\ f ds ¥z{s)r{s,t)Oy\ 
zcx Jr 



Now a Picard iteration for the related quantity 

\KT{r,t)Oy\\ 



C x {r) := sup 



I /CI 



(21) 



(22) 
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Figure 2: In the proof of Theorem 1, we need to bound a 
sum over all paths of length n starting in X, the support 
of tCx, and ending in Y, the support of Oy', see Eq. (24). 
A path corresponds to a sequence of local Liouville terms 
(izt) with overlapping supports. In the depicted situation of 
a two-dimensional lattice with nearest-neighbor interaction, 
path C\ would contribute to the sum for n = 10 and C2 for 
n — 5. To simplify the calculation for upper bounds, the sums 
are extended to contain all paths starting in Y (if Vol(F) < 
Vol(X); all paths starting from X, otherwise). Hence, in the 
bound for n = 5, also paths like C3 are taken into account. 



can be used to obtain a bound for ||G(r)||. Inserting 
Eq. (21) in Eq. (22) gives 

C x {r) < C x {t) + 2 sup \\t z {s)\\ f dsC z (s), 
C x (t)<S(X,Y)\\0 Y \\, (23) 



where S(X, Y) = 1 for X n Y ^ and S(X, Y) = 0, 
otherwise. The second line follows from JCxOy = for 
Liouvillians K x € L x with X n Y = 0, and ||/COy|| < 
||/C|| 1 1 Of 1 1 in general. Starting the Picard iteration for 
C x (r) with Eq. (23) and Z := X leads to 



Cx(r)<\\0 Y \\J2 



(t - r) n 



with 



n=0 



III 



(24) 



c «= E E ■•• E S(Z n ,Y)Y[ sup JlzM 

Z 1 CZ Z 2 CZ 1 Z n CZ„-! 



8 = 1 



se[r,t] 



Now we can exploit that the Liouville terms are of finite 
range a and that they induce a finite maximum num- 
ber Z = ma,xz-.e z ^o Vol(Z) of nearest neighbors; Sect. II. 
The sum in Eq. (24) runs over all paths from X to Y. 
Depending on whether Vol(X) or Yo\(Y) is larger, the 
number of such paths with length n can be bounded by 
the number of all length-n paths starting in X or 7, re- 
spectively See Fig. 2. This gives the simple bound 



c n < V x ,Y(2\e\) n , and thus, 

00 An 

C x (t)<Vx.Y\\0 Y \\ J2 



(25) 



where 6 := (t - r)Z\£\, D = \d(X,Y)/a], and ZV x ,y 
is the minimum of the numbers of Liouville terms izi 
supported in X and Y, V x .y = min{Vol(X), Vo\(Y)}/Z. 
We have also used that c n = for all n < D, as one 
needs at least D Liouville terms of overlapping support 
to pass from the subsystem X to subsystem Y, such that 
5(Z n ,Y) ^ 0. Using induction, the sum in Eq. (25) can 



be bounded by Y£Ld TH — e 
Theorem 1 follows, 

||G(t)||<Vx,y||/Cx||||0, 



-D. 



see Appx. C. Hence, 



,(t-r)Z\t\e-D 



VI. CONCLUSION 

We have shown that the evolution of an observable 
with support Y under a quantum master equation with a 
short-range Liouvillian can be approximated by the evo- 
lution with respect to the truncation of the Liouvillian to 
a subsystem V D Y. The error decreases exponentially 
in the distance of Y from the complement of V. With 
this tool, we derived an error bound for Trotter decom- 
positions of the propagator. Those results correspond to 
efficient simulation techniques for open-system dynamics 
on classical and quantum computers and provide rigorous 
bounds to finite-size effects. 



Acknowledgments 

We gratefully acknowledge inspiring discussions with 
J. Eisert, C. Gogolin, V. Nesme, and T. J. Osborne. This 
work has been supported by the EU (Minos, Qessence), 
the EURYI, the BMBF (QuOReP), and the Studiens- 
tiftung des Deutschen Volkes. 



Appendix A: Operator and super-operator norms 

In this work, two of the Schatten p-norms [42] are em- 
ployed. The oo-norm of an operator O G B(H) is defined 
as its largest singular value and is equal to the operator 
norm, 



\0\\ 



inn- -Ml 

\0\\ ■= sup 

Wen WWW 



(Al) 



where |||V')| = \J (ip\ip) denotes the vector 2-norm. The 
oo-norm is the physically relevant norm for observables. 
The 1-norm, of an operator E B(H) is defined as the 
sum of its singular values and is equal to the trace norm, 



lOlli = HOIItr ^TrV0W. 



(A2) 



It is the physically relevant norm for states, i.e., den- 
sity matrices [43]. Those operator norms induce corre- 
sponding norms for super-operators T G B(B(H)). The 
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(oo^oo) -norm is defined as 

IIS! := WTW^ := sup 
and the (1— >T) -norm is 
IIT 



\TO\\ 
oeB{u) ||0|| o 



ll-s-l • 



Ijrojji 

= sup . 
OeB(-H) ll^lll 



(A3) 



(A4) 



In order to switch between the Schrodinger and the 
Heisenberg picture, one needs to consider the adjoint 
of a super-operator T, defined by 



(A,TB) KS = (T^A,B) HS y A . BeBm , 



(A5) 



where (•, -)hs denotes the Hilbcrt-Schmidt inner product 
(A,B) HS := Tr(A^B). The (1 l)-norm is dual to the 
(oo^oo)-norm in the sense that 



imi c 



= sup WTOW^ 

l|0|loc=l 

= sup sup \(X, TO)hs| 

lloiU=i ||x|| 1= i 

= sup sup \(T^X,0)ns\ 
ll-xriu =1 lion =i 



= sup 

11*11,=! 



\T*X\ 



lrt| 



i->i 



(A6) 



This allows us to relate the appropriate norm of a prop- 
agator T in the Heisenberg picture to the norm of the 
corresponding propagator in the Schrodinger picture. 
For more on properties of the norms, see, for example, 
Refs. [42, 44]. 



Appendix B: Properties of the propagators 

The derivative of a propagator with respect to its sec- 
ond time argument is given by 



d t T V (s,t) = T V (s,t)C v (t). 



(Bl) 



Using the defining properties d s Ty(s,t) = 
—Cv{s)Tv{s,t) and rv(t,t) = id, Eq. (Bl) follows 
from the equation 

= 8 t id = d t [T V (t,s)T V (s,t)] 

= [d t T V (t,s)]T V (s,t) +T V (t, s)[d t T V (s,t)] 

after applying Ty{s,t) to it. 

Let us explain why the propagators are norm- 
decreasing, i.e., 

Ms,t)0\\ < \\0\\ \/£eh A ,s<t,OeB(H). (B2) 

The adjoint propagator r^(s,t) (see Appx. A), describes 
the time-evolution in the Schrodinger picture, p(t) — 
(s,t)p(s), where p(t) denotes the system state at time 



t. First of all, we note that r^(s,i) is a completely posi- 
tive, trace-preserving (CPT) map since it can be written 
as a product integral [45] , 

rt(*,t)= lim TTe^)^. 

3 

Every factor e c ^ tjS)Atj is an exponential of a constant 
Liouvillian and is hence CPT. Thus, the finite products 
are CPT maps and, since the set of CPT maps is closed, 
also the limit T\s,t) is a CPT map. Then Eq. (B2) fol- 
lows from the norm duality ||T|| = WT^^ = \\T%^ 
[Eq. (A6)] and ||Tt||i_»i = 1 for all CPT maps T+. The 
latter has, for example, been shown in Ref . [27] . 



Appendix C: Bound on the partial exponential sum 



In the following, we prove that 



Eoo 
n— 



N n \ 



< e 



xe-N 



\fx > 0, N € N . 



(CI) 



=:/jv(aO 



Note first that, for N = 0, 

f (x) = e x < e xe = g (x) V x > . 
The functions f N and g N obey the differential equations 

d x fN+i{x) = fN(x), d x g N+1 (x) = g N (x) V Xi jv- 
For all N > 0, the initial values /jv(0) = and ajv(0) 



-N 



obviously obey /jv(0) < <7iv(0) Vjv>o- Consequently, 
/iv(a;) < 9n(x) V x >o implies fN+i{x) < g N +i{x) V x > . 
This proves Eq. (CI) inductively. 



Appendix D: Bound on a sum of exponentials 

In the following, it is shown that 

oo 

53 n K e- n < 2eD K e-° Vk>0,D> 2k+1 E N. (Dl) 

n=D 

Due to the definition T(a, x) := f 00 dtt a ~ 1 e~ t of the in- 
complete Gamma function, one has 



n K e- n <r(/s + l,D-l). 



(D2) 



n=D 



The bound 

T(a, x) < Bx a - X e- X Va > 1, B > 1, x > 
of Natalini and Palumbo [46] , reads for the choice B = 2 

r(o, x) < 2x a - 1 e- x Va > 1, x > 2(a - 1). 
Together with Eq. (D2) one obtains 

oo 

53 n K e- n < 2(D - l) K e- D+1 V K>0 , d>2 k+ i 

n=D 

and hence the acclaimed Eq. (Dl). 
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Appendix E: Trotter expansion of a propagator 

For two times q < t, we derive the Trotter error bound 
\\Mq,t)- J] rz{q,t)\\Y 

ZcV 

< (t - q) 2 ZYo\{V)\£\ 2 e^ q)W (El) 

that is employed in the proof of Theorem 3. To this 
purpose, let us determine an upper bound for the right- 
hand side of 

\\rpq,t rjiq,trjiq.t\\ ^ \\rpq.t rpq,trpq,t\\ 

\\ 1 C+£~ 1 C 1 £ lly - \\ 1 C+£ ~ 1 C 1 i II' 

where denotes the propagator for a Liouvillian 

IC(t) e L, C(t) e L obeys the preconditions of Theo- 
rem 3, and £(t) e hz is a local Liouvillian term with 
support Z. We denote the inverse of a propagator T^ 1 



by T^. Using d q T q / = -K{q)T%\ d t T q / = T q /lC(t), 
T^T^ = T£ , T£ — id, and applying the fundamental 
theorem of calculus twice, one finds 

rpq,t rpq.trpq.t frpq.t rpt.qrpt.q • i\rpq,trpq,t 

1 C+l^ 1 C 1 £ ~\ 1 C+l 1 l 1 C ~ m ) 1 C 1 l 

= f ds d s (T^T^T^TfTf 

Jq 

= f dsT q ^ e (C(s) -T°' q £( S )T^ s )T°' q T^T^ 

Jq 

= f ds f drT^d r (T s t ' r C{s)Tl> s )T*' q T s /Tf 

J q J q 

= fds / S drr^,T;-[^(r),£( s )]r;-^ t r^ t . 

J q J q 



The time arguments occurring in the integrand are or- 
dered according to q < r < s < t. The norm of the 



A bound for the norm of 



propagators is \\T£ || = 1 M s < t - A Dounc 
the inverse propagators can be obtained from their repre- 
sentations as time-ordered exponentials [27, 45] , yielding 

^ ex P(X* dr l|£(*)ll) V ^<t- With those properties, 
the triangle inequality, and the norm submultiplicativity, 



\rpq,t rpq.trpq,tU < - 

\ 1 c+i ~~ 1 c 1 e II — 



ds / dr \\[£(r),C(s)\\ > 

q Jq 



(s-q)\l\ 



< (t-q) 2 Z\e\ 2 e«-«M. 



This bound and the inequality ||7iT 2 - TiT 2 || < 
ll^i II 11^2 - T 2 \\ + ||Ti - Ti||||T 2 || can now be used itcra- 
tively, separating one local propagator Tf'* after another. 
As Cy is a sum of Vol(V) terms £z, Eq. (El) follows. 
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